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Abstract 

! A theory for studying the dynamic scahng properties of branes and relativistic topological defect 

o 

\ networks is presented. The theory, based on a relativistic version of the level set method, well- 

J>-,' known in other contexts, possesses self-similar "scaling" solutions, for which one can calculate many 

quantities of interest. Here, the length and area densities of cosmic strings and domain walls are 
\ calculated in Minkowski space, and radiation, matter, and curvature-dominated FRW cosmologies 

with 2 and 3 space dimensions. The scaling exponents agree the naive ones based on dimensional 
. analysis, except for cosmic strings in 3-dimensional Minkowski space, which are predicted to have 

(N 
t — 

f-*) I densities are a factor of approximately 2 lower than results from numerical simulations of classical 

(N 

■ field theories. An expression for the length density of strings in the condensed matter literature is 



a logarithmic correction to the naive scaling form. The scaling amplitudes of the length and area 



corrected. 
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I. INTRODUCTION 



The solutions to some of the most interesting problems in physics depend on a better 
understanding of the dynamics of fields far from thermal equilibrium, particularly in particle 
cosmology, where we seek mechanisms for generating a baryon asymmetry for generating 
density fluctuations P], and perhaps for generating primordial magnetic fields P]. Signif- 
icant advances have been made recently in studying non-equilibrium dynamics of phase 
transitions, both in theoretically (see e.g. |^] for a review) and numerically, where we can 
now perform real time simulations of a quench with leading thermal corrections included 



|. One aspect is still not yet well understood: the approach to equilibrium after a phase 
transitions of a field theories with topological defects. 

At the same time, the last few years has seen an explosion in theories involving various 
kinds of extended objects or branes, both solitonic (like topological defects in field theory) 
and fundamental . Most of the interest has lain in special configurations of branes of various 
dimensions, and the spectrum of states in those backgrounds. However, an interesting new 
scenario has emerged in which the Universe began with the branes in thermal equilibrium, 
the brane gas Universe 0] 

Both branes and topological defects in relativistic field theories obey the same equation 
of motion (at least for configurations with curvature small compared to the inverse width 
or fundamental scale), and so it clear that both may be discussed at the same time. Hence 
the theory presented in this paper can be applied to both brane gases and networks of 
topological defects. The general technique is independent of the space-time dimension and 
the codimension of the brane, but quantitative predictions must be taken case by case. The 
cases worked out in detail here concerned defects of codimension 1 and 2 in Friedmann- 
Robert son- Walker (FRW) spacetimes of dimension d = 3 and 4. 

It is believed that when extended topological defects are formed, self-similar or scaling 
behaviour emerges at large times, in which a characteristic length scale of the field configu- 
ration, ^, increases with time as a power law: 

Dynamic scaling can be seen in the order parameter of many condensed matter systems 
undergoing rapid quenches, and there are now quite sophisticated techniques for calculating 
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correlation functions of the order parameter 



. They fall into two classes. Firstly, there 



are those based on a large expansion, where N is the number of components of the order 
parameter, which are applicable to Ginzburg-Landau theories. The second is applicable to 
systems with extended topological defects, in which the order parameter (p obeys an equation 
equation oc 6F[(f)\/6(j), where F is the Ginzburg-Landau free energy. Allen and Cahn |9| 
proposed that the velocity of defects marking a phase boundary was proportional their local 
mean curvature. This proposal, now termed motion by mean curvature, was later rigorously 
proved |lo| . 

Relativistic scalar field theories with spontaneously broken global symmetries (Goldstone 
models) also exhibit dynamic scaling. Significant progress has been made on the theory 

nnn n 

of O(A^) scalar field theories at large A^, both classical [IJ, |15|, |16[ and quantum [Ij] (at 
large A^ the leading order in the quantum theory is the same as the classical theory). These 
works have established a theoretical basis for the scaling observed in numerical simulations 
HQQ. The theo.. has a.o bee. used to ea.e„.ate ...owave a„d density 

fluctions. To date, however, analytic approaches to the dynamics of topological defects are 
few. 



There are several numerical simulations which broadly support the dynamic sca^ 



pothesis for topological defects, including domain walls 



and global strings 



27 



21 



22 



23|, gauge strings 



ing hv 



25,^ 



I . All the simulations are consistent with the linear scaling law over 
the range of the simulations, although Press, Ryden and Spergel suggested that the results 
for domain walls would be better fitted by ^ ~ t/ln(t); however, more recent simulations 
with a larger dynamic range [2^ are not consistent with the logarithm. 

There are also string simulations based on direct integration of the equations of motion of 
one-dimensional objects, obeying the Nambu-Goto equations, which may be derived as the 
first approximation in an expansion in powers of string worldsheet curvature 

HI 0,0 They 

do not include any way for energy to be lost from the network, but if one considers "infinite" 
strings only (strings longer than the causal horizon size), an approximately linear scaling 



law is found 1^ 1^ • However, the simulations are plagued by kinkiness persisting at the 
resolution of the simulation, associated with the production of small loops of string, which 
does not appear to scale. It has been suggested that this is because the natural length scale 
for loop production is in fact the string width, where loops would become indistinguishable 
from large amplitude oscillations in the field 2^- Indeed, numerical simulations of the fields 
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appear to support this hypothesis, although the latter authors suggest that the 
"protoloops" in their simulation are in fact a transient effect. 

A programme to understand analytically the results of the Nambu-Goto simulations has 



developed over the years 



33, 



37 



38, 



In its simplest form, the model parametrises 



the string with one length scale ^, which is defined from the invariant length density of 
infinite string Lqo through Lqo = 1/^^- This length density can change in two ways: through 
stretching as the strings participate in the Hubble expansion, and through loop production. 
Loop production is parametrised by the so-called chopping efficiency c, the fraction of string 
lost to the network in the timescale ^. The Hubble stretching depends on the mean square 
string velocity v^. The phenomenological equation is then 



i = Hil + v')^ + c/2. 



(1) 



Further work introduced two other length scales to describe the correlation length and 
the inter-kink distance. However, there are many unknown parameters in the model which 
greatly restricts its predictive power, despite attempts to measure them 3^. A different 
approach was adopted by Martins and Shellard 39| who promoted the r.m.s. string velocity 
f to a time-dependent parameter, to model the reduced rate of loop production of slower 
strings. The velocity-dependent one-scale model equations are (neglecting frictional terms) 



-2Hv + k{l ~ v^)C\ 



(2) 



where c and k are, in the simplest version, constants. It is this velocity- dependent one-scale 
model which 2^ use to make their claim that the production of loops on the scale of the 
string width seen in field theory simulations is a transient. 

In this paper a potentially far more powerful analytic technique for describing the motion 



of strings is developed. The technique was outlined in |40|, |4l| and applied to relativistic 
domain walls in 2 and 3 space dimensions. It is here further extended into a partial treat- 
ment of p-branes in D space dimensions, and fully applied to relativistic strings in 3 space 
dimensions. It is based on the w-theory of Ohta, Jasnow and Kawasaki (OJK) 1^, and its 
descendents |431], which describes the motion of defects obeying the Allen-Cahn equation. 
The relativistic generalisation of the Allen-Cahn equation is the Nambu-Goto equation, 
in which, loosely speaking, the acceleration of the defect is proportional to its local cur- 
vature, with proportionality constant c^, where c is the speed of light. More precisely, the 



Nambu-Goto equation is equivalent to the requirement that the world- volume of the p-brane 



embedded in the d-dimensional spacetime has zero extrinsic curvature. How a 
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ose 
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V defects 



4J. The 



derived from a field theory obey this equation is a matter for debate [30! 
theoretical approach develops systematic expansions of the geometrical equations obeyed by 
the defect world-volumes in powers of the width divided by the local curvature, which reduce 
to the Nambu-Goto equation in the limit of small curvature. The approach of Arodz 3^ 
makes it particularly clear that the Nambu-Goto equation is really a consistency condition 
for a smooth defect-like solution to exist. 

It is therefore plausible that we can forget about the details of the field theory and concen- 
trate instead on the properties of extremal (zero extrinsic curvature) surfaces embedded in 
higher dimensions. If one finds such surfaces, then provided their curvature is small enough 
one can be confident that there is a solution of the field equations representing a smooth de- 
fect centred on that surface. A formalism for studying extremal, and more general, surfaces 
has been developed over the years by Carter [4^, which makes clear the geometrical nature 
of the Nambu-Goto equations through close attention to the tensorial properties. 

The present approach we introduce scalar fields with the intention that the loci of 
constant should be extremal surfaces: these are the level sets of the title. The fields can 
also be interpreted as coordinates normal to the brane surface: in this sense the approach 
can be thought of as orthogonal to Carter's. We derive the equations that the must 
satisfy, which are non-linear, and so therefore do not seem to represent an improvement on 
the original field theory or the Nambu-Goto equations. However, one can derive equations 
for surfaces which are on average extremal, when we average the fields with a Gaussian 
probability distribution. With this Gaussian ansatz, one can also calculate analytically 
important quantities, such as the brane or defect density. 

The results for {D — l)-branes (domain walls) are extremely encouraging when compared 



to the numerical simulations 



21 



47L l48| . The theory predicts a scaling law for the area 
density in 3 dimensions, but not only does it predict the scaling exponent, it also predicts 
the scaling amplitude to within a factor of about 2, which is not bad given the approximations 
made. The prediction for {D — 2)-branes in 3 dimensions (strings) is also challenging: the 
theory gives a logarithmic scaling violation in Minkowski space, with the length density 
deoending on conformal time 1] as log{ri)/ri'^. Looking for such scaling violations will be a 
good way to test the theory, although computationally very challenging. 
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The theory also describes the behaviour of defects formed from initial conditions with a 

nnn 

slight bias in the expectation value of the field favouring one vacuum over another |241I471I48|. 



It is found that the defects disappear exponentially fast at a critical conformal time r]c, which 
scales with the initial bias f/ as r/c ~ U^^^. Indeed, part of the motivation for this work was 
o account for this kind of behaviour observed in simulations by Coulson, Lalak and Ovrut 



47[ and Larsson, Sarkar and White 



3- 

Finally, in making comparisons with similar results in the condensed matter literature, 
an expression for the length density of strings in 3 space dimensions in the condensed matter 
is corrected (see Section IV Dj) . 

In this paper we shall work a conformally flat (i- dimensional Friedmann- Robertson- Walker 
space-time with coordinates , such that d = D + 1. The metric is given by 

gpiv = a^(r7)diag(-l, 6^j), (3) 

where t] is conformal time, giving an affine connection 

TP,. = {5',6l + 6PX-g,.g^'){a/a). (4) 



II. FIELD EQUATIONS 

In this section we shall firstly study model field equations for topological defects of codi- 
mension = 1 and N = 2, which correspond to walls and strings respectively in D = 3. 
We shall see that we can find approximate solutions to the field equations near surfaces of 
codimension N which have zero extrinsic curvature, and whose other curvature radii are 



large compared with the widt 



shown in various ways in 



viatn oi tn 

3, 31, 3 



" the defect. These results are well known and have been 



but the approach here is slightly different and worth 



exhibiting in some detail for the later sections of the paper. 



A. Domain walls 



Let us first consider a theory with a single scalar field (p, with action 



S 



J d''xV^i^^d,<pdy + v{<p)y (5) 
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from which we derive the field equation 



^ ._d,{^g^-^d.)<P+^^0 (6) 



We shall suppose that the potential V has the symmetry —(p, and moreover that its 
minima are at = ±v, with V{±v) = 0. If we impose the boundary conditions: 

— oo) = —V, (f){x^ — >• +oo) = +v, (7) 

and make the ansatz 

d^<P{x)=0 {i^ = 0,...,D-l), (8) 

then the theory has a one parameter family of domain wall solutions, with = at a;-^ = X^. 
If the potential is quartic, 

v{<p) = li<i>'-vr, (9) 

then the solutions are 

= 't;tanh[M(x-° - X^)], (10) 

where M = \^v. Thus the width of the defect is controlled by the parameter M^^. The 
defect can be thought of as centred at X^, where the field vanishes, with a width parameter 



B. Strings 

The simplest theory to exhibit string-like solutions is the Abelian Higgs model, which 
has action 

S = -J d'^xV^ [If^uF'^-' + D^<p*Dy + , (11) 

where is a complex scalar field with covariant derivative = 9^0— ieA^0. The potential 
V is taken to respect a U(l) symmetry cf) — > e*"0, with a circle of minima at |0| = v. If we 
impose the boundary conditions in the {x^^^,x^} plane 

0(r ^ oo) = ve'\ (12) 

where = {x^~^y + {x^^ and tan^ = x^/x^~^, then by continuity must vanish some- 
where in the plane. If we furthermore assume translational invariance in the other d — 2 
directions in spacetime, we find a two-parameter family of static string solutions, labelled by 
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the coordinates of the centre of the string, {X^ ^,X^}. In the radial gauge Ar = these 
solutions take the form 

4>{x) = /(p)e*'^, Ai = —ipia{myp), A^ = 0, (13) 

ep 

where pi = (x — X)^~^, P2 = {x — X)^, p^ = (pi)^ + (p2)^, tan(^ = P2/P1, and if'' is the 
unit azimuthal vector in the {x^~^,x^} plane. These solutions cannot generally be found 
analytically, even when the potential has the renormalisable and gauge invariant form 

v{^) = lxm'-vy2r. (14) 

However, they are easily found numerically, and exhibit similar properties to the domain 
wall in that away from the centre of the defect the fields approach their vacuum values 
exponentially, at rates controlled by the masses of the fields = ^/Xv and m„ = ev. 
Defining a dimensionless coordinate z = m^p, and (3 = [nis/m^Y = A/e^, one has jl2] 

/ ~ 1 - hz-^'^ exp(-y^2), a ~ 1 - aiz^'"^ exp(-2). (15) 

In the case /3 > 4, the asymptotic form of / is 1 — z~^ exp(— 22;). 

Again, the string can be thought of as centred at {X^~^,X^}, with thickness m^, al- 
though for light scalars (/3 <^ 1) there is a thicker scalar core where the scalar field asymptotes 
to its vacuum value. 



C. Solutions in curvilinear coordinates 



These are however rather special solutions with a high degree of symmetry. Let us instead 
look for (if necessary approximate) solutions, corresponding to defects centred on a more 
general surface X^(cr"), with a = . . . ,p = D — N . We choose a new set of coordinates 

= {a°',u^}, where A = 1, . . . , N, with the intention that the equations of the surfaces 
can be written 

u^{x) = 0. (16) 
We write the metric in these new coordinates 

^ _ I daX- dpX daX ■ ObX \ _( lo-fi ^ Ba 

y Oax ■ dpx Oax ■ ObX j y Na/s Gab 
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where the dot indicates a contraction with respect to the original metric g^^p. We may choose 
the coordinates so that, at least at u"^ = 0, the and a" are locally orthogonal, or 

^A/3L.=o = 0- (18) 

In fact, with walls and strings in D = 3, these are only 3 or 4 conditions on the metric 
respectively, so we know we can make a coordinate transformation so that this is true 
everywhere, and not just at u"^ = 0. 

Note that the upper left (p + 1) x (p + 1) block of G^y, denoted 7q,/3 in Eq. (fT7|). is the 
embedding metric on surfaces of constant u^, which they acquire by virtue of being surfaces 
embedded in a spacetime with metric g^^. 

We can also write the inverse metric 



We define 



(19) 



n 

^ jJLU 



[Nap hAB + NApNB^ 



(21) 



I^AB ^ Q^A . Q^B^ (20) 

and use the convention that the indices a, /3, etc. are raised and lowered with 7^/3 and 7"^ 
(defined as the matrix inverse), and that the indices A, B, etc. are raised and lowered with 
h^^ and its matrix inverse Hab- Hence 

jBa 

One can show that 

det G^''' = det 7"^ det h^^ , (22) 

and hence that G = ■y/h, where G = det G^^, 7 = det 7^/3, and h = det h^^. 

We have two projectors associated with the constant surfaces, one which projects onto 
the surface and the other which projects onto the subspace spanned by the vectors d^u^. 

P|f, = y'^^d^X^dpXp, Pi, = hABd^u^dpu'', (23) 

withP|jt + P^, = 5^ 

Let us study the field equation for the theory of an A^-component scalar field $ in these 
new coordinates: 

^ dJV^y^^dp)^ = ^dA(^/^h^^dB]^-'%-. (24) 
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u(x) < 




u(x) > 



FIG. 1: The coordinates = {u", n"^}, where a = 0, . . . ,p and A = 1, . . . ,N, which are chosen 
so that = will be the extremal surface on which the topological defect sits. Illustrated is a 
1-brane in 2 + 1 space-time dimensions, located at X^(a) = x^{a, 0). 



where Oa = Oa — Na^Oo (where we use notation mirroring that of Moss and Shiiki 
At the surface m"^ = it is possible to make a coordinate transformation amongst the u"^ 
coordinates so that they are orthonormal, that is, h^^ = 6^^. 

This choice of coordinates is different from the one used in other works on solving defect 
equations of motion in curvilinear coordinates where coordinates jo"", p } are 

constructed away from the surface by setting 

x'^(a", p^) = X^(a") + pM(^^°), (25) 



where = Oax^I^a^q, and ua ■ ub = Sab- The coordinates coincide only when A^^ 



Carter 



0. 



45| also uses orthonormal vectors in the surface, and is careful to express quantities as 



space-time tensors. Table H] contains a summary which compare his notation and conventions 
with this work. 

In contrast to previous work, it is here more convenient to use the unnormalised d^u^ as 
basis vectors, as we are interested in the surfaces generated by gaussian random fields m^, 
with unconstrained derivatives at = 0. 

We now try to find approximate solutions to Eq. ()24|) . A promising avenue is to look for 
solutions which are independent of a", in which case Eq. fj24j) becomes 

where is the extrinsic curvature of the constant hyper surf aces, given by 

1 



(26) 



-7 



-7). 



(27) 
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This work 


Name 


Carter 


Relationship 


PL 


First fundamental tensor 




PL = 




Orthogonal projector 




PC — 1 M 


daX^ 

(a = 0, . . . ,p) 


Tangent vector 


(A = 0,...,p) 


iA^Lftn = riAR. 


(A = l,...,iV) 


Normal vector 


{X = l,...,N) 


Ax^Ay^ = 6xY-, 




Second fundamental tensor 




— K Pd 



TABLE I: Comparison of notation and conventions with that of Carter 



45|. 



The ansatz $ = ^{u^) can only be self-consistent if both K"^ and h^^ are independent of 
cr". This is still a difficult equation to solve, so the next step is to look near surfaces where 
the extrinsic curvature vanishes. Transforming to the orthonormal coordinates (|^^ near 
those surfaces we have the approximate equations 



d d ^ dV ^ 
$ -I ~ 



(28) 



By "near" we mean the region where If^^^Syi^l ^ Eq. |2H1 is solved by ^{p^), 

the original defect profile. Hence we are guaranteed approximate solutions to the field 
equations near smooth = (extremal) surfaces. The argument in this section can be 
straightforwardly extended to gauge fields and so the task of solving the field equations has 
been replaced by the task of finding extremal surfaces. 

The extrinsic curvature will generically vanish only at = 0, and be non-zero 
elsewhere in spacetime, and so the static solutions $ will not be exact. However, we should 
be able to find approximate solutions $ = $ + yj, with the perturbation being sourced by 
the departures of from zero: 



(29) 



where ^r^S,,^') is the retarded Green's function for the scalar field fluctuation operator, 
given by 

[-Wn + V"(0)]A^(e,r) = m - a, (30) 

with Aji{^, = for C,^ < If the extrinsic curvature decreases with time, the source for 
the perturbation ip dies away, and we should not have to worry that our initial assumption 
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that \ip\ ^ |$| is rendered invalid. In fact, the dynamic scahng hypothesis holds that 
K"^ ~ where ^ is the average curvature radius of the defect network. 

There are in fact special cases for which K"^ = everywhere, and exact curved defect 
solutions exist. These represent travelling waves on walls and strings js^, although they do 
not obey a superposition principle because of the non-linearities in the field theory. 

This brings us close to the controversial subject of radiation from defect networks. We 
postpone this discussion until Section IVTl 



III. EXTREMAL SURFACES 

We saw in the last section that if we could find a suitable surface of constant (which 
without loss of generality we can choose to be = 0) satisfying = 0, an approximate 
solution of the field equations could be found. We shall now derive the equations that 
must satisfy in order that u^{X) = be an extremal surface. 

Differentiating once with respect to the world- volume coordinates a", we find 

df,X^d^u^{X) = 0. (31) 

(This equation is of course true independently of the choice of the coordinates Using the 
embedding metric we can covariantly differentiate (j31|) by acting with (—7)~^^^9q,(— 7)^/^7°"^, 
where 7 = det'ja/s, to obtain 

(p+i) ox^d^u"^ + ^"^d^X^dpX-'d^d^u^ = 0. (32) 

The operator 

(P+l)n = (_^)-V25^(_^)l/2^a/35^ (33) 

is the covariant d'Alembertian in the surface = 0. 

The equations of motion are obtained by extremizing the invariant area of the surface 



11 



12 



Anv[X] = J <F+\^-^{X) (34) 
with respect to the enbedding coordinates X^{a). The result is 

(p+i) nx/^ + r(;^7"^5aX"9/3X^ = 0, (35) 

where is the affine connection derived from the metric g^j,^. 
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The reader will notice the appearance of the tangential projector P|j in equations 
which we replace by g^'^ — P^ Combining and to eleminate the d'Alembertian 
we find 



g^^^ - hABd^u^d''u''\ {d.d.u'^ - TP^^d.u^) = 0. (36) 

This is the fundamental equation of motion for the fields u^{x), which strictly only applies 
at = 0. 

The equations also follow from a variational procedure. Using the fact that G = 'y/h, 
and that G = g, one can show that the invariant area of a p-brane can be re-expressed in 
terms of the as 

Ai^^[u^] = J d^x^Vh5^{u). (37) 
Varying with respect to and dividing by y^— 7 gives us 



-5''{u)h^V^, ih^-^hABd^u'') = 0. (38) 



This can be shown to be equivalent to 

S''{u)PfV,d,u'' = 0, (39) 

and hence Eq. (j36|) at = 0. In orthonormal coordinates, for which h-^^ = 6ab, Eq- f!38p 
becomes 

{u)V^n^^' = 5 ^ (m) i^^ = , (40) 

where is the extrinsic curvature. Thus we can identify K^^ = P^^V adyU"^ as the extrinsic 
curvature tensor, or equivalently the second fundamental tensor (see Table E] and 145^). 

The restriction that the equations apply only at m"^ = complicates the finding of solu- 
tions, and we assume that we can extend the equation K"^ = to all u^. It is not obvious 
that non-trivial solutions exist to the extended equations, because such a solution would 
be a foliation of space-time in which all leaves have zero extrinsic curvature. As mentioned 
above, some non-trivial solutions are known 571 but there is no general existence proof 

Q.J 
. However, we could equally well look for solutions to 

= f{u^), with f{u) any function which vanishes at u = 0, so there should be a certain 

amount of freedom. Furthermore, we will be looking only for perturbative solutions to the 

extended equations. 
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IV. AVERAGE EXTREMAL SURFACES 



The equations of motion ()36p are not easy to solve, as they are non-hnear. However, 
they have the distinct advantages over the alternatives. The equations of motion for the 
coordinates of the u"^ = surfaces (jS^ are non-local: defects generically self-intersect. 
This non-locality generally defeats analytic approaches, and also makes numerical simu- 
lations algorithmically difficult, as one must devise an efficient scheme for searching for 



self-intersections ISc 



34 



35l |. The equations of motion for the underlying field theory are 



also non-linear, and in the gauge of the Abelian Higgs model (and other gauge theories) 
have a gauge covariance, which precludes the naive application of techniques like large N. 
Numerical simulations of field theories are relatively straightforward, but require significant 
amounts memory to allow the scale of the network to grow much larger than the width of 
the defect. 

Instead of trying to find families of surfaces whose curvature is exactly zero, we shall 
find surfaces who curvature is zero on average. The average will be taken with respect to a 
Gaussian probability distribution for u^. We assume that the distribution function remains 
Gaussian throughout the evolution, which is similar to the approximation underlying the 
large approximation in scalar field theory. Indeed, we should expect there to be a similar 
large limit in this theory. 



A. Gaussian averaging 

Our starting point is an ensemble of coordinate functions u^{x) with an assumed Gaussian 
distribution. Thus the average value of all observables of interest Q{u^, d^u^), which we take 
to be functions of and its derivative d^u"^, are evaluated with probabihty distribution 

dP[u^] = Pwexp (-i j d^xd''yu\x)C2U^,y)u''{y)^ , (41) 

where C^^{x, y) is the 2-point correlation function. 

We are often interested in densities, which means that the observable Q is evaluated at 
a particular point x. This means we can simplify the evaluation of the averages from a 
functional integral to an ordinary one, as we now demonstrate. 
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First, let us take the Fourier transform of the observable, 

n{u\x\ = I fc)e^'-^«+^^^^^"^(^). (42) 

We now introduce current densities La{x) and K'^{x), according to 

La{x) = IaS\x-x), K'^ix) =k'^S'^{x-x), (43) 
so that the expectation value of Q{x) is given by 

{n{u\dy)) = J J rfPMf2(/,A;)e^/'^'*^(^-(^)-^-^-(^))""(^). (44) 

Performing the integral of the random field u^, we find 

(fi("^aX)> = / (0 (^s^C.*)^-*^-^."'-"'^'"*''''-*"'". (45) 

Substituting the form of the functions La and Ka from (PHj) . we find 
(27r)^ (27r) 

where 



C^^iv) = limC^^(x,i/), d.C^'^iv) = hm-^C^^{x,y), 

dAC--i,) = limA|_c- (.,,), (47) 

and 7] = x^ . We note that we expect correlation functions to be spatially homogeneous 
but to depend on conformal time non-trivially, reflecting the symmetries of the background 
spacetime: hence the explicit conformal time dependence of the two-point correlators eval- 
uated at the same two points. At this point we recall that the Fourier transform of the 
observable VL may be written 

n{l,k)= J d^'ud^^'^TT l](n^, n^)e-''^^^-''>^ . (48) 

We can economise slightly on the length of expressions by introdcuing some new notation. 
Defining N{d + 1) -dimensional objects j and / by 

3 = {lA,k^A}, / = K,<}, (49) 
with a scalar product (j, /) = IaU'^ + ^a^t^! we can write 

m)) = I d^^^'^'^f (1^^ e-i(-^^)-(-/), (50) 
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where C is the covariance matrix 

Finally, we may perform the integrations over the Fourier transform variables j to obtain 

= [(27r)-^('^+^)/2][detC]-5 J ^^(^^+1)/ n(f) e-5(/.c-V). (52) 
Hence the average of the observable can be found with an ordinary integral, as claimed. 

B. The covariance matrix C 

In our case the covariance matrix C is restricted by the assumed FRW form of the back- 
ground. It will be homogeneous and isotropic, but not time-independent. We will also 
assume an 0(iV) symmetry between the N coordinate functions u^{x). Hence, the basic 
two-point correlation function at zero separation may be written 

C^^(r;) = = 5^''C(r;), (53) 

We shall also define a function M^,, from the two-point correlator of d^u^: 

d,d.C^^{r]) = {dy{x)d.u^{x)) = 5^^M,,{r^). (54) 

The assumed spatial isotropy of the distribution function dictates the form of M^i,: 

( T{r}) 

\^ SmnS{ri) 

With this definition it is not hard to show that S{ri) = —C"{ri), where C"{rj) = 

Two-point correlators with odd numbers of derivatives also occur, as the ensemble is not 
time translation invariant. The correlator with one derivative is 

d,C^^{ri) ^ {dy{x)u^{x)) = \5^^8lC{r^l (56) 

and with three, 

{dy{x)d,dpu^{x)) = i,.p{r))5''^. (57) 



(55) 
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Again, symmetry restricts the form of 'y^up'- 



7000(^?) = 2^(^)' lOmniv) 



1 



(58) 



It is interesting to note that ^i^^p = \{.M^^^p + M^p^^ - M^p^^). 
Thus the covariance matrix can be written 



f C \ 

he T 



6rr,,nS 



(59) 



Its inverse is easily found, and defining the determinant of the upper 2x2 block A 
(TC — ^C*^), we can write 



c- 



JAB 



( 



\ 



T 



he 



\ 



he c 







1 

A 



(60) 



^mn g j 



The determinant factor in the probability distribution is also straightforward: 

[detC]-^ = [5^A]-^/l (g^) 

Often, we will want to find expectation values which are independent of SqM"^, mainly because 
the integrals are easier to evaluate. By integrating over -Kq one can easily show that 

(n(«^,a,n^)) = 



[(27r)^5^C]-^/2 / ci^Mci^^Tr fi(u^, t:^) e"^ 



(62) 



It is very convenient to rescale the integration variables in the probability distribution, 

A 



u'^ — » u^^/C and nf TTf\fS, in which case 



(n(M^,a;,M^)) = 



(63) 



C. Averaging the null extrinsic curvature condition 

The averaging procedure is greatly aided by rewriting the equations of motion (j36|) in the 
following form: 



-9 



d 



-g det h 



pu 



0. 



(64) 
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The procedure now is to linearise the equations of motion by taking the Gaussian average, 
and then to find a self-consistent solution for the fields u'^{x, rj). We will require the following 
identities, which are proved in Appendix lEl 

(det h d^d.u^) = (det h)d^d,u^ + ^7^^, [9^^^^^ ^-^^ (^^) 
(det h dpu"") = (det h)d,u'' + ^ (^^(det h)^ M^,9,n^. (66) 

The expectation value of the determinants in ()(i5|(ifi|l can be expressed in terms of the two- 
point correlator M^^ (defined in Eq. 1^ 

(det h) = m n M,^,^f^^^5^^i::^^-5:,\::z, (67) 

where 5pl"'p^ is the identity tensor in the space of rank N antisymmetric tensors, defined 
in Appendix^ The right hand side of (jU7j) resembles a determinant, and we introduce the 
notation detM to refer to it. We can also define a kind of cofactor for M^j,, 

W = m g''^''^6^p':Y,n7^::::M,,,, . . . M.^^jd^tM. (68) 

Putting the pieces together we find that the linearised equations for surfaces which are on 
average extremal are 

2 



d_ 

dgp" 



{g^" - g^'g^'^T^ ) detM [d.Au'' + -m"%.^.9.m^ 

2^, 



- [dru'' + -M"'Ma.9am^)] = (69) 
With the assumed symmetries for the correlation functions, these equations have the form 

where ^{rj) and v depend on T, S*, and the background cosmology parametrized by a, and 
must be taken on a case-by-case basis for each A^. 

D. Linearised equations for walls and strings 

In our three-dimensional Universe, the cases of most interest are A^ = 1 (domain walls) 
and N = 2 (gauge strings). A^ = 3 corresponds to gauge monopoles, which do not scale 
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For N = 1, (det h) = M^^g^"", while for AT = 2 (det h) = iKM^^g^"")^ - M^^M^"]. 
We then find, for FRW backgrounds (see Appendices ID ![ lEll ) 

^(^) = I -^S/S) + a{v) [D - 3 (T/S)] {N = 1), 

^ 1 -(2/(/^ - l)MS/S) + aiv) [{D - 1) - 4 (T/S)] (iV = 2), 

where 0(77) = rja/a, and 

[ [D-2-2(T/S)]/D (iV = 2). 

In scahng solutions, we expect S and T to have power law behaviour, and so as long as we 
are not near a transition in the equation of state of the Universe (such as that between the 
radiation- and matter-dominated eras), fi and are constant. Thus, imposing the boundary 
condition that m"" be regular as 77 ^ 0, (fTOj) has the simple solution 

where A'^ 2T(z/ + l)-u^(?7i) as — 0, and (1 — /i)^/4 = z/^. The form of the initial power 
spectrum is taken to be a power law, with index q, and an upper cut-off at \k\ = A. 

We may now evaluate T/S and f^, and self-consistently solve for the undetermined pa- 
rameter fi. It turns out that one must take u = —(1 — /u)/2 if all the integrals are required 
not to diverge as A ^ 00. This also gives regular solutions as r/ — ^, because as it turns out, 
(x > I. With this choice, C scales as r/~(-^+'J)^ g g^^fj ^~{D+q+2)^ 

In the following, we will take the power spectrum to be white noise, g = 0, as is consistent 
with a causal origin for the defects in a phase transition. There the power spectrum of the 
scalar field from which the defects are made has a g = power spectrum at long wavelengths, 
and so we should take the fields to have a similar power spectrum if we want to reproduce 
the statistics of the defects from the statistics of the zeros of u^. 

Using standard integrals of Bessel functions, and defining the parameter f3 = 2u — D — l = 
H — D — 2, we find (see Appendix in|l 

1 ^ J 2(D+2+/3) /y^N 

S I (g+2)(g-2) /AT-O^ ^ ^ 
L 3(D+2)+2/3 ~ 

provided /3 > 0, which ensures that the integrals for 5* and T are defined. Given the 
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N 


a 


/5 


(TIS) 




1 





p + 2) 


D-1 
4 


3 D-1 

4 4 




oo 


(D+2)(D-3) 
2D 


D 

3 


2D-3 
3D 


2 





(D+2)(3-D) 
D-1 


(D-2)(D-1) 
D+3 


(D-2)(D+5) 
D(D+3) 




oo 





D-2 
3 


D-2 
3D 



TABLE II: Values for parameters /3, and T jS of the self-consistent solution to the linearised 
equations of motion ()7U() for the N fields n"^ for = 1 (domain walls). In the special cases of 
Minkowksi space (a = 0), and curvature-dominated FRW cosmologies (a = oo), exact values can 
be found for all D. 



expressions for T j S \s, easy to show that 

f (D-l)(D+2+2/3 f/V = n 
^2 ^ I 2D(D+2+/3) '^"^ /^g^j 

I (D-2)(D+2+2/3) /AT -9^1 ^ 
L D[3(D+2)+2/3] ~ 

To find /3, we must solve the equations derived from (f7T|): 

^^|«[^-3(T/5)] + (I^ + 2) (Ar = l), ^^^^ 

[ a - 2 - 3(T/5)] + (D + 2)(3 - D)/{D - 1) (iV = 2), 

which are quickly seen to be quadratic. One can obtain results in simple closed form in 
Minkowski space {a = 0) and curvature dominated universes {a = oo) which are displayed 
in Table |n] For other backgrounds the solutions may be written down in closed form, but 
are not particularly illuminating as they are fairly lengthy expressions. 

Instead, numerical values of (3, T/S and v"^ for particular cases of interest are given: 
radiation-dominated (a = 1) and matter-dominated (a = 2) 2 and 3-dimensional universes 
(Tables Uni 113) • 

Note that for strings (A^ = 2) in 3 dimensions in Minkowski space (a = 0), for which 
fi = D + 2, and hence (3 = 0, which does not satisfy the requirement (3 > for the integrals 
defining 3(7]) and T{ri) to be convergent. One finds that a logarithmic scaling violation 
appears, and S,T (x log{Ar])r]~^^~^'^~^'^\ We also have a solution with /5 = 0, and therefore 
logarithmically divergent S and T, for walls in 3-dimensional curvature-dominated universes. 
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L> = 3 


D = 2 


a 


/9 


(TIS) 




P 


(T/S) 







5 


1/2 


1/2 


4 


1/4 


3/8 


1 


6.72 


0.43 


0.52 


5.36 


0.21 


0.39 


2 


8.83 


0.36 


0.55 


6.90 


0.18 


0.41 


CXD 





1 


1/3 


-1 


2/3 


1/6 



TABLE III: Values for parameters /5, and T/S of the self-consistent solution to the linearised 
equations of motion (|7()j) for the N fields n"^ for = 1 (domain walls) in D = 2, 3. Values listed 
are for Minkowksi space (a = 0), radiation-dominated (q = 1), matter-dominated (q = 2), and 
curvature-dominated FRW cosmologies {a = oo). 





D = 3 


a 


{T/S) 




P 





1/3 


1/9 





1 


0.22 


0.14 


3.65 


2 


0.20 


0.20 


4.75 


oo 


1/3 


1/9 






TABLE IV: Values for parameters /?, and T/S of the self-consistent solution to the linearised 
equations of motion ()70() for the N fields ti"^ for A = 2 in D = 3 (strings). Values listed are for 
Minkowksi space (a = 0), radiation-dominated (a = 1), matter-dominated (a = 2), and curvature- 
dominated FRW cosmologies (a = co). 

V. AREA DENSITIES FOR WALLS AND STRINGS 

Armed with the mean-field solution for u^[x) we can now calculate anything that can 
be expressed in terms of local functions of the field and its derivatives, provided of course 
that we are able to perform the Gaussian integrals involved. Here we derive formulae for 
the area densities of defects, where by "area" we mean the world-volume of the {p + 1)- 
dimensional hypersurface u"^ = 0, which has dimensions of (Length)"^. We must be careful 
to distinguish between various kinds of area: there is invariant or proper area which is a 
coordinate-independent quantity, and there is also the projected p-dimensional area. The 



21 



latter quantity is what one would obtain by simply measuring the p-dimensional area of 
the defects at a particular time. This quantity is the most convenient to calculate for 
comparison with numerical simulations, which is a good thing as the proper area density 
is far harder to calculate. One must also bear in mind that area densities are coordinate- 
dependent quantities: in the cosmological setting we will need to convert between comoving 
area density and physical area density by multiplying by the appropriate power of the scale 
factor a, which is a~^. 

Here we give figures for the projected area densities of walls and strings in D = 3. 
They can be compared with results from numerical simulations of the field theories and give 
surprisingly good agreement given the uncontrolled nature of the approximations made. 

A. Proper area density 

The proper area density ^ of a p-dimensional defect in D space dimensions is 

AUx) = J d^^'a' ^5\x - X{a'))/^g. (77) 

Making the coordinate transformation from to E,^ = {a°',u^} near the world volume of 
the defect, we have 

AUO = I d^^'a' V^S^^\a - a')5''iu^)/^^G, (78) 
Recalling the results of Section pi C^ . we can perform the integration over a', to obtain 

^P^ = 5^(M^)|det/i^^|i/2^ (79) 
where the reader is reminded that 

= d.^u'^d^u^g^" (80) 

Thus the problem of calculating the proper area density is reduced to finding the Gaussian 
average of A^^ in ()79p. The conversion factor from comoving to physical area is given as 

with N = D -p. 
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B. Projected area density 

Easier to measure and to calculate is the projected area density, which is defined as 

A:'d = I <P^' v^5^(a; - X{a'))/^, (82) 

where is the spatial part of the metric. The induced D-dimensional metric on the 
p-dimensional surface m"^ = is 

iDab = daX^d^X^gm,, (83) 
where a,b = 1, . . . ,p. As for the proper area density, one can show that 

A^ = 5^(M^)|det/l^^|l/^ (84) 

where 

h""^ = diU^u^'gt (85) 

Note that g^j!) is defined as the matrix inverse of goij, and is not the spatial part of g^^'^. The 
conversion between physical and comoving area is again 

A'n^^ys = a-^'Al. (86) 



C. Average projected area density: walls 

We can now use the averaging formula ()63p to find the mean value of the operator A, 
which when we specialise to domain walls (A^ = 1) gives 

{A^-') = (^/l / dud''^ e-^"^-^-^'^-^ (87) 

The integrals are easily performed to give 



^-y^ T{D/2) ' ^^^^ 



4^. This is the 



a well-known result originally derived by Ohta, Jasnow and Kawasaki 
comowmg' projected area density: to obtain the physical projected area density, one multiplies 
by a^^: 
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D. Average projected area density: strings 



For strings (A^ = 2), the average we need to calculate is 

where the rescaled quantity h^^ is given by 

h^^ = Trfrrf. (90) 

Now, 

det/i^^ = le^cesD/^^^/^^'', (91) 



-(^AC^BBTcfTcfTcfn^, (92) 



which suggests that we construct the following antisymmetric matrix: 



= nfiifeAB, (93) 

such that 

det/i^^ = ^^^-. (94) 

Thus in order to calculate the average area, we need the probability distribution for fij. 
At this point we specialise to D = 3, as the calculations are considerably simplified by 
introducing the vector 

4>k = 2^ijkfij, (95) 

whereupon 

det/i^^ = (96) 

The probability distribution for = |0fc0fc|^ is derived in Appendix |F[ and turns out to be 
remarkably simple, giving 

{AD = dW<j> 5^(«^)0e-^"""^-^e-^. (97) 

A simple calculation now shows that the comoving projected length density for strings in 
D = 3 is 

(^3) = (98) 
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Note that this disagrees with the formula derived by Toyoki and Honda |49| , but agrees with 



Scherrer and Vilenkin 



Toyoki and Honda write the 3D string length density as 



A\ = 6{u^)6{u^)\Vu^ X Vm^I = 6{u^)6{u^)\Vu^\\Vu'^\ cosOu (99) 

where 612 is the angle between the vectors Vu^ and Vm^. They then average 612 over a 
uniform distribution, separately from and u"^, which is incorrect. 



E. Projected area density: higher 



Scherrer and Vilenkin |55[ used an elegant argument to derive their value for the projected 
area densities of walls, strings and monopoles in D = 3, which can be generalised to any 
and D. They noted that a string was located at the intersection of two surfaces = and 
u"^ = 0, and therefore the length density string could found by computing the length per 
unit area of the lines of = in the surface = 0, and then multiplying by the area per 
unit volume the surface = 0. That is, 

Al = AlAl (100) 

which clearly has the correct dimensions. One can easily check that this gives the correct 
result Al = (S/ttC). It is immediately obvious how to generalise the formula to any D and 
N: 

n=p 



Thus 



where N = D — p, 



Av (^\"'" np + i)/2] .^Q2) 

^^"Ucj T[{D-N + l)/2\ ^^^^^ 



F. Quantitative results 

It is shown in Appendix ((GJ that 

S_ _ 1 D + 2 + I3P + 1 
C ~ ^ 4^^2 

In the special cases of = 1, 2 one can substitute for from Eq. (f73j) to obtain 



(103) 



^ I 1 D(D+2+/3)(l+/3) ^/v-n 

c I DmD+2)+2m+p) (N-'?\ 

4(D-2)/3 
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D= 


-3 — ^ \/^ 


D 


_9 41 _ 1 /S^ 

z y 


a 


Theory 


Simulation 


Theory 


Simulation 





1.91r/-i 


0.88(0.14).r/-i-00{0-03) 


l.llr?-! 


0.77(0.23).r/-0-99{0-03) 


1 


2.02r/-i 


0.93(0.13).r/"0-99{0-0i) 




0.93(0.17).r/-i-00{0.02) 


2 


2.1677-1 


0.96(0.12).r7-i-00{o.oi) 


1.24?7-i 


1.15(0.23).r7-0-99(o.oi) 



TABLE V: Comparison between theoretical and numerical simulation values of the domain wall 
defect scaling density in Minkowski space, FRW radiation and FRW matter dominated universes 
(a = 0, 1, 2 respectively) in 2 and 3 dimensions. The numerical values are taken from ,29]. 



It was shown in Section IV El that the projected area density is proportional to (S/C)^^"^, 
and therefore classical scaling behaviour for all defects is predicted, unless (3 = 0. By 
classical scaling, we mean that the area density goes in proportion to conformal time as 
naive dimensional analysis would predict: a p-dimensional area density in D dimensions 
should be proportional to f]^'^ , as indeed it is in this theory. When /3 = 0, as is the 
case for {D — 2)-branes in D = 3 (strings) in Minkowski space, and for {D — l)-branes in 
curvature-dominated FRW backgrounds, logarithmic violations to naive scaling appear. 

We are also able to compute the scaling amplitudes, the coefficients of the relations 
between the area density and the appropriate power of time. These can then be compared 
with numerical simulations. The scaling projected comoving area densities for walls and 
strings, in the radiation and matter eras are displayed in Tables |3 IVll Note that in Table 

a the results for strings in matter and radiation-dominated universes have ben taken from 
, who give proper area densities. These have been converted to projected area densities 
by dividing by ((1 — f^)-^/^), where v is the average speed of the string. While not strictly 
the correct procedure, it gives a good enough answer given the uncertainty. 

To convert between comoving and physical areas, one uses the formula Aphys(t) = 
a~^A{ri), and the fact that a{Tj)rj = (1 + a)t. 

The scaling amplitudes differ from those obtained in numerical simulations of 0^ theory 



j29| and of the Abelian Higgs model 



25, 2a, by a factor of about 2. However, it should be 
noted that there are large errors on the central value. The authors of Refs. 

HQ did not 

look for logarithmic scaling violations in the area density for strings in Minkowski space, 
choosing instead to fit to a simple power law. Finding such a violation is numerically very 
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D=3 A\ = 


_ 5 
~ ttC 


a 


Theory 


Simulation 





3.6r/-2 iog(7;A) 


(11 ± 1)t/-2 


1 


6.8r?-2 


(18±6)t?-2 


2 




(14 ± 4)??-2 



TABLE VI: Comparison between theoretical and numerical simulation values of the string scaling 
density in Minkowski space, FRW radiation and FRW matter dominated universes (a = 0,1,2 



25| and 



26l |. with the latter 



respectively) in 3 dimensions. The numerical values are taken from 
converted from proper to projected area densities. The numerical fits in Minkowski space did not 
look for a logarithmic scaling violation. 

demandine;, large dynamic range is required. 



Biased initial conditions 



One may also ask how the network behaves when a small bias is introduced into the 
initial conditions, that is, if (u^(xi)) = U^. In numerical experiments simulating biased 



initial conditions for strings [51[ it is found that as the bias is increased the string passes 
through a transition from a phase with a finite fraction of percolating "infinite" string 
and with a power-law size distribution of loops, to one without infinite string, and with an 



exponential size distribution for the loops. In numerical simulations of domain walls 
it is found that even for very small initial biases, for which the walls percolate, the system 
still evolves away from the percolating state and eventually the large walls break up and 
disappear. Similar behaviour is well-known in in the study of quenches of condensed matter 



systems with a non-conserved order parameter [42 , 

The theoretical description of this behaviour is fairly straightforward. Introducing a bias 
into the initial conditions for walls alters the Gaussian average of Section IV CI to 

1 [J 



and hence 



(27r)'^/2 



dud^TT S{u)\Tri \ e 



■Uu-U/VCf-^7T,5*i7Tj 



I S T[iD + l)/2] 
TiC T{D/2) 
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(105) 



(106) 



It is clear that this form is common to all defects in all dimensions: if {A]^)o is the unbiased 
average area density, then the result of including a bias is 

{A'o) = {A],)oe-'^^'/^, (107) 

with an obvious generalisation to > 1. If the system is close to being self-similar at some 
initial time r]i when the magnitude of the bias is U and the fluctuation around that value 
C{rii), then one can predict that the area density goes as 

A ~ r/-^/2 exp(-cf/2r7^), (108) 

where c is a constant. One can also show that the time rj^ at which the defect density falls 
to a fraction of its scaling value as 

Vc = Vi{uV2C{ri,)y'^'' . (109) 

The simulations by Larsson and White are consistent with ()l()(i|l and (I109|l in D = 2, but 
do not have sufficiently good statistics in D = 3 \^\. Coulson et al. [47] did not attempt a 
fit of the form p06p to their simulations. 



VI. SCALING AND ENERGY LOSS 



There is an apparent inconsistency in our conclusions for topological defect networks. We 
started by establishing that one could find approximate solutions to the field equations by 
finding extremal surfaces in spacetime, and then constructing static solutions in coordinates 
which moved with the surface. We then showed that one could construct random surfaces 
in FRW spacetimes which are on average extremal, whose average area density obeyed a 
classical scaling law with conformal time 1]. The assumption is that there are defect-like 
solutions which are somehow close to static solutions centred on these random surfaces. 

There is a problem with this picture: the defect area density decreases with time and 
therefore the energy in the form of defects also decreases. This energy must go somewhere, 
and an obvious channel is into propagating modes of the fields, or radiation. However, 
it is difficult to reconcile the idea that the network energy is lost into radiation with the 
perturbative approach to finding curved defect solutions, which assumes that the deviation 
from the comoving static solution decreases with the curvature of the defect. 
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Indeed, there is good numerical evidence that the perturbative approach works in certain 



cases 



44| . The configurations where it has been tested are coUiding travelhng waves, 



either sinusoidal 



26j | or more complex J^- When travelling waves are correctly prepared 
to the recipe laid down by Vachaspati 57J, the collision does produce perturbations in the 
form of radiation, which is however exponentially suppressed with decreasing curvature. 

It should be noted however, that pure travelling waves are obtained from very special 
initial conditions. A random defect network is not prepared so carefully and it appears that 



it does radiate by an as yet poorly understood mechanism |25l. l2q. The radiation shows no 
sign of being exponentially suppressed with increasing curvature. What is clear is that one 
or more of the assumptions implicit in the perturbative approach to finding curved defect 
solutions must be violated. Two possibilities are that the extrinsic curvature is much larger 
than maybe due to kinks, or there are non-linear radiative processes, perhaps involving 
the breather modes jssl . 

VII. SUMMARY AND CONCLUSIONS 

To summarise, this paper describes a new analytic technique for describing the dynamics 
of a random network of branes or topological defects, applicable to the brane gas universe 
or a cosmological phase transition. It is a relativistic version of a well-known approach 
in condensed matter physics, due to Ohta, Jasnow and Kawasaki [4^, which uses a mean 
field approach to find approximate solutions to the AUen-Cahn equation for the motion of 
a surface representing a phase boundary. In the relativistic version, the surfaces are branes 
or defects obeying the Nambu-Goto equation (i.e. they have zero extrinsic curvature), but 
the condensed matter analogues can be obtained as a certain limit (see Appendices ID H and 
IE which acts as a check. In rederiving these condensed matter results an expression for 
the length density of strings due to Toyoki and Honda j3| has been corrected (see Section 



EM- 

In most cases the prediction is that the (generalised) area density of a p-dimensional defect 
in D dimensions should scale with conformal time as r/~*^^~P\ with a scaling amplitude of 
0(1). This appears to agree quantitatively with numerical simulations of domain walls 
|29l |48[. In certain cases, such as strings in D = 3, there is a prediction of a logarithmic 
violation of the naive scaling law. There are further predictions for defects with biased initial 
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conditions, for strings in 3D, and for {D — 1)- and {D — 2)-branes which would be interesting 
to test. 

From the point of view of the brane gas Universe, it would be interesting to look at 1-, 
2- and 5- branes in higher dimensions. One of the most interesting features of the brane gas 
scenario is that it offers and explanation of why the Universe has three large dimensions: 
strings do not generically interact with each other in more than three dimensions, and so 
winding modes can never decay. It is only a 3-dimensional subspace, where the winding 
modes can interact with each other and annihilate, which can expand and becaome large. 
It follows from this idea that strings cannot scale in more than 3 space dimensions, as 
there is no opportunity for the initial winding modes to break up into closed loops in the 
conventional picture of energy loss by a string network. It is therefore important to see 
whether the theoretical techniques presented in this paper predict scaling for strings in 
higher dimensions. 
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APPENDIX A: PRELIMINARIES 



Define a projector onto the rank N antisymmetric tensors (which is also an identity 
operator for those tensors). 

= {^Z • • • C + ''Sned perms on u,) (Al) 
This projector has the properties 

_ }^ i/i...z/jvM]v+i---/id ( A9^ 

"/xi...Atjv jyi^^ _ jY'jpMi---AiJvMJV+i---Md'= ' K^^J 

S:ui..MN Xpl-.pN — Xpi-PN ( A'i\ 

r/' 

v-Miv m{d-N)\ ^ ' 
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Define the matrix h"^^ = d^u^dyU^ g^'^ . Then 

det h = ^^ g^^-^ . . . g^---d,,u^^d,,u^^ . . . d.^u''^ d,^u^^ eA,...A^eB,...B^ (A5) 
Define the antisymmetric rank tensor 

i^w-Miv = ^m"^' • • • d^^u^'^eA^.^Ar,- (A6) 

Then we may write 



1 

Note that 



det/i = — F^,.„.,F^— (A7) 



= d^u'^^dyu^' ■ hA,B^ det h. (A8) 



APPENDIX B: AVERAGING THE EXTREMAL SURFACE EQUATION 

In a general spacetime, the equation for a D-dimensional surface with zero extrinsic 
curvature is 

{g^" - hABd^u^d'u'') {d^d^u^^ - vi.dru'^) = o, (Bi) 

where 

hAB = (A-'"')-' 

The surfaces of constant m*" satisfying this equation have = 0. Note that the following 
is the projector onto the tangent space of the surface of constant u'^: 

= g^^^^ _ hAsd^u^d^u^ (B2) 

Thus, if we write v'" = du^ as the coordinate vectors normal to the surfaces of constant m*", 
we can express the equation as 

Pf^f^Vu = 0. (B3) 
Recalling the identity |X8l we see that the following equation holds: 



det h (d^dyu^ - Tl^drU^)] = (B4) 



Hence, in order to obtain the equations for surfaces whose average extrinsic curvature is 
zero, we need to average the quantities ^iQihd^dvVp and det/iS^w'". 
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1. The Gaussian average {deihdfj,di,vP) 

Exploiting its antisymmetry, we may rewrite the tensor as 

F,,...,, = N\5^::Zd.,v} . . . (b5) 

Hence the determinant becomes 

det h = NW^\::2,d^''v}d,y . . . d'''' d^^u^ . (B6) 

We introduce m^^ = d^u^dnU^ (with no imphed summation), which is an unnormahsed 
projector orthogonal to the surfaces of constant u"^. Then 

det h = NW^-^-^---ml^,^ . . . m;r,,^. (B7) 

Hence 

{deihd^d^u^) = {det h)d^d^u^ 

with no implied summation on the index C, and the wide hat symbol is used to denote a 
term removed from the product inside the angle brackets. We now use the relations 

{Oyd^u'') = M^,(^^«, (B8) 

{dpu%d,u^) = 7pM.5^'', (B9) 

("^Jii.i^M^^'^^) = Ifiini^d^iU^ + -f^^^^d^^u'^' (BIO) 

from which we can immediately derive 

(det h) = iVl^'^-'^-'^i-'^-M^,,, . . . M^,,, (Bll) 

and 

^ ^ {deth)^p^,d,u^ (B12) 



NdMp, 



2. The Gaussian average (det hdru'~^) 
It follows from the previous section that 

+2m5^'^-^-^'-^^M^,,, . . . aQ,^ . . . M^^^^M^^rd.^u''. (B13) 
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3. More definitions 

We define a kind of determinant det tlirougli tlie relation 

d^M = NW-^'^'''-''^M^,,, . . . M^^,^ = (det h). (B14) 
We can tlierefore define a cofactor for M^,y, wliicli we denote M^'^, tlirougli 

W = NW^'^-^'^'"'^-''^M^,,, . . . M^^.J^m (B15) 

Thus we may write 

(det h d^d.u^) = d^M (d^d.u^ + ^W-ip^.d^u'^^ (B16) 



and 



(det h dru^) = detM ( d^d^u^ + ^M^^Mprd^u^ ] (B17) 



APPENDIX C: THE MEAN FIELD ZERO CURVATURE EQUATION 

Putting the results of Appendix together, we find that the Gaussian averaged equa- 
tions for zero extrinsic curvature surfaces is 

^g'-'' - d^M [dpd.u'' + ^M'^"7«^.5ru^ 



(CI) 



For future convenience we will break this equation down into four terms: 

Ta = (g"'' - ^] d^Mdf^d^u^ (C2) 



dg, 



Tb = ^ [g^"" - ^) detMM^^7,^,9AM^' (C3) 
Tc = (g'-' - d^MTld.u'' (C4) 



^g^^,. 



Td = ^ (^g'" - ^) detMM^^M«,ri,9.w^ (C5) 

Before reducing this equation further in cases of definite we will need the following explicit 
expressions for the correlation functions M^i, and 'j^fiu, consistent with the spatial 0(D) 



33 



symmetry: 

(to, 

M,, = \ (C6) 

V S6:^ 

7000 = -^T, 7omn = —-SSmn, ImOn = 1m.nO = -SSmn (^7) 

Note that 

l^up = \ (M^,,p + M^p,, - M,p,p) . (C8) 

We also need the Christoffel symbol for a flat FRW background, which has metric g^^y = 
o?'{vi)rj^y, where rj^y = diag(— 1, 1, 1, 1) is the Minkowski metric. It is 



rj. = ^ + sX - (C9) 



APPENDIX D: ZERO CURVATURE FOR N = \ 

The simplest case is with one coordinate field m, which is appropriate for domain walls 
in 3 dimensions. Here we have 

detM = gP^Mp, = M (Dl) 
detMM"^ = g""^ (D2) 
detMM^^M^^ = M^r (D3) 

The required derivatives with respect to the metric are also easily found: 

d 

-detM = M>"' (D4) 



dgf,u 

detMM''^ = g^'^'g''^ (D5) 



dg^v 

Using the exphcit form of M^y we can also write down 

M = a-'^{-T + DS) (D6) 

a. Term Ta 

Using (jnH) and (jDH) we find 

TA = {g^''M-M^'')d^dyU. (D7) 
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Using ()C6|) and ()D6|) this simplifies to 

Ta = -a-^DS 

b. Term Tb 

Using (|D2|l and (|D4|l we have that 



Using ()C7|) one can quickly show that 



Tb = 2a~^DSu. 



c. Term Tq 

Using (jnH), (jnH) and one finds 



= [M(1-D)^°^-2M°^ + Mc/°^]9am 
^DS[D - 2 - (T/S)]m. 



d. Term Td 

Using (jnil) and (jDS)) we find that 



'a 



= 2[(2 - D)g''g^^ - g^'g'^ - g^'g'^]M,x 
= 2[(1 - D)M°" - ^°"M] 

Substituting the known forms of M^'^ and M we arrive at 

Td = 2a-^DS[l - {T/S)]u. 
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1. The equation for u 



We can now construct the Gaussian averaged or "mean field" equations of motion satisfied 
by the coordinate function u, which is apphcable to domain walls when D — 3. The equations 
are made from the four terms we calculated in the previous section: Ta + Tb — Tc — Td — 0. 
Putting them all together, and dividing by the factor a~'^DS, we find 

ri 

with 

„^[D- ?,{TIS)] (r^^) - 2 (^r^l) , (D15) 
= [{D - 1) - {T/S)]/D. (D16) 

We can recover the well-known AUen-Cahn equation for the overdamped motion of domain 
walls by identifying the damping constant F = a/d, and neglecting T/S, riS/S, and the 
second order time derivative of u: 

ii = T^^^V'^u. (D17) 
APPENDIX E: ZERO CURVATURE FOR = 2 

When N = 2 the expressions for the various quantities involving M in the equations of 
motion are still straightforward to evaluate: 

ditM = {g'^'^'g'^'''^ - g'''''^g'^'''')M^,,,M^,,, 

= {M^ - M'^'M^^) (El) 



detMM"^ = {g^'^g^^"^ - g'"'^g^'^^)Mf 



= g'^^M-M'^^ (E2) 
ditMM^'M,, = {g-^g^^^^ - g^^-g^-^)M^^,,M,, 

= M^M - M'^^M^r: (E3) 
We also need to differentiate two of these expressions with respect to the metric g^i,. 

detM = 2(M'''^M - Mf'^Mx") (E4) 



dgi^u 

^ detMM"^ = {g^'^'g^^M + g'^^M^"' - g^'^'M"^ - g^^M^""). (E5) 



dg^^u 
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Introducing a further piece of notation, that M ■ M = M^^M^ we can show 

M -M = a-^{T^ + DS^) 
detM = M^ -M -M = a-^DS\D - 1) - 2T / S)] 

Idt/s 



1 



[{D - 1) - {T/S)]5„ 



a. Term Ta 
Using lE4l we find 



Ta= (M^-M- M)d'^u^ - 2{M'"'M - M^'^Mx'')d^d^ 
Hence, using (jETj) and lE8|l . 

Ta = a-^D{D - 1)S^ 



b. Term Tb 

Using dill) and we find firstly that 



Using (jlSl) and (jEg]) we find 



ditMM"\,^,9AM^ = -a-'^DSKD - 2)5 - T]m^. 



Putting the two expressions together we find 



c. Term Tq 

From ()E1|) and ()E4|) we can immediately write down 

Tc = [g'^'iM^ -M-M)- 2M'"'M + 2M'"'M^'') T^^^dxu 
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c 



Using jCgl) we find 

Tc = a"2[(l - D){M^ - M ■ M) + AM°^{M^° - M) + 2(M^ - M ■ M)] f ^ ) . (E16) 



With equations ()E1|) and ()E8|) we arrive at 



Tc = -a-'^D{D - l)S\D - 3) - 2{T / S)] Pj m^. (E17) 



(i. Term To 

For this last term we begin with 



9 ^ 

Hence from (IE3I) we see that 



A'l'pA _ Q _ 



Z^) (^) r . (E18) 



(7'^^detMM""M,Ar;:,9AU^ = a'^g^'^M - M°°)Mo°(l - D) if 

= a-^D(D - 1)ST (^^) M^. (E19) 

The second term in expression D is more comphcated. From ()E5j) and ()(^9|1 we have 

- g^^'M''^ - g''^M^'')M^yj:l^. (E20) 

After some algebra we find 



detMM^^M^ArJ^ = [^"^(M^ - M ■ M) - 2M°^M + 2M^aM°^] 



= a~^D{D - l)S^5l.) (E21) 
Subtracting ()E2Hl multiplied by dru'" from ()E2H) . we arrive at 

Td = -a-''D{D - l)S^[l - {T/S)] if. (E22) 

1. The equation for {N = 2) 

We can now construct the Gaussian averaged or "mean field" equations of motion satisfied 
by the coordinate functions in the case N = 2, appropriatre for strings in 3 spatial 
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dimensions. The equations are made from the four terms we calculated in the previous 
section: Ta + Tb — Tc — To — 0. Putting them all together, and dividing by the common 
factor a-^D{D - 1)S\ we find 

V 

with 

H^[D-2- 3{T/S)] (ry^) - (^rj^^ , (E24) 
= [{D - 2) - 2{T/S)\/D. (E25) 

We can recover the results of Toyoki and Honda for the motion of overdamped strings in 
D = 3 by setting their diffusion constant F = a/ a, and neglecting T/S and rjS/S. In this 
case we get 

uC ^ £v2«^, (E26) 
3 

which is identical to their equation (3.10). 

APPENDIX F: PROBABILITY DISTRIBUTION FOR Fij 

The definition of the antisymmetric tensor F^j is 

Fij = diU^u^'eAB. (Fl) 

The probability ditribution for Fij is therefore constructed from the Gaussian probability 
distribution of diU^. Fij is antisymmetric, so we need only consider half of the non-zero 
elements, e.g. by imposing i < j. Moreover, it is convenient to scale out the variance of 
diW^, defining variables irf and fij as follows: 

diU^^y/STif, F,^ = Sf,, (F2) 

where 

{diu'^{x)d,u^{x)) = S{t)5ijS^^. (F3) 
Hence, the probability distribution for is 

Pih)\i<j = / n ^^-'^^-^ih - ^MeAB)U<r (F4) 
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Using the Fourier representation of the 5-function, 

P(m,^, = / (gp /nig^e-'^-.'-E..'-<'«-'^'f-', (F5) 

where P = D{D — 1) /2 \s the dimension of kij. 

We now do the ttJ^ integrations in turn, starting with the highest A. First, note that 

^(/«)i.<, = / / n ge-i-."'.'--'.'-?, (F6) 

where there is now no restriction on the sum over i,j in the second exponentiaL Second, 
define the variable — k'^^nj. Then we have to evaluate the integral 



A 



(2vr) 



Doing the tt/ integral first, this is 



I{h,) = I ^e-^'^^^^^ (F8) 



(27r) 

= det-5M (F9) 

where 

Mij = 6ij + kikkjk- (FIO) 
At this point we specialise to 3D, where we can write 

hj = ^ijkPk- (Fll) 

Hence 

M,j = 6,,{l+p')-p,pj. (F12) 
The eigenvalues of this matrix are 1 (twice) and 1, so 

det-^M = {l+p^)-\ (F13) 

Thus the probability distribution of fij is 

In 3D we can replace fij by (pk = Hi<j^^^^ fiji and the integral may be easily evaluated to 
give 



where 0^ = 0^0^. 
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APPENDIX G: INTEGRAL FORMULAE AND CORRELATION FUNCTIONS 



In this appendix we perform the integrations necessary to evaluate the functions C, S 
and T, defined in Section, which we repeat here for convenience. 

5^^C(r/) = (n^(a;)n^(a;)), 



with 

/ 



(Gl) 



\^ 6mnS{r]) 

We shall also evaluate C for the mixed correlator (d^u'^u^). We recall from Eq. ()73|) that 
the linearised solution for with the correct boundary conditions is 

with z/ = ±(1 — /i)/2. If we demand regularity and convergent integrals as ?7 — 0, we must 
take the negative sign here, as it will turn out that yU > 1. 

In order to calculate the two-point functions it is useful to define the following integral 

POO 

I{p,a,T)= dzz'PJ^{z)Jr{z), (G3) 

^0 

which has the value 



1 r(p)r r+v+\ 

/(p, a, t) = -, , ^ -V, (G4) 

' ^ p p~cr+T+l \ Y ( p+cr+T+l \ Y ( p+o-~T+l \ ' ^ 



provided Re(cr + r + 1) > Re(p) > 0. The first inequality comes from the condition that 
the integral be defined as z ^ 0, and the second from requiring that it converge as z oo. 
There is a simple pole in at p = 0. We can see that this comes from the z^^/^ behaviour of 
the Bessel functions as 2; — > oo, and corresponds to a logarithmically divergent integral. 
Defining the Fourier transform of the correlator C in the usual way through 

Ciri) = J £^Ck{ri), (G5) 



we see from the solutions for that 

1 

[vr])^ (27r 



C(V) = zAhTt!^ / dzz^~'-^^Jliz)PAk) (G6) 
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where z = kvrj, and fin = 2tt^/^/T{D/2) is the volume element of a (D — l)-sphere. We 
assume a power- law form for the power spectrum of A^, 

where A is a high wavenumber cut-off, satisfying Avrj ^ 1 for all rj of interest, and ai is the 
variance. Hence, defining i3 = 2u — D — 1 — q, 



Let us now calculate S from 



DSiv) = I ^k'Cj^iv)- (G9) 



One can straightforwardly show that 



DS= ^[^^ , — -/(/?, z/,z/). (GIO) 



The correlation function T is obtained from 



Given the identity 



one can show that 



d_ ( Ju{z) \ _ Ju+i{z) 
dz\ z'' ~ z" 



(G12) 



T=—^^^——^—^ — -Up u + + (G13) 

Note that the ratios S/C and T/S depend on the initial conditions only though the power 
q, which appears in (3; 

S 1 u, v) 



C D{vriy I{2 + p,u,u) 



(G14) 



S I{/3, V 

A little more algebra shows that 

S_ _ 1 £' + 2 + /?/3 + l 



(G16) 
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Note that the ratio S/C appears to have a simple pole at /3 = 0: however, when the cut-off 
is in place this is replaced by a logarithm, with 

^^^\og{Kv^). (G18) 
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